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A survey is given of recent developments on the resummed small- a; evolution, in a frame- 
work based on the renormalization group equation, of non-singlet and singlet structure 
functions in both unpolarized and polarized deep-inelastic scattering. The available re- 
summed anomalous dimensions are discussed for all these cases, and the most important 
analytic and numerical results are compiled. The quantitative effects of these small-x 
resummations on the evolution of the various parton densities and structure functions 
are presented, and their present uncertainties are investigated. An application to QED 
radiative corrections is given. 
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A survey is given of recent developments on the resummed small- a; evolution, in a framework based on the 
renormalization group equation, of non-singlet and singlet structure functions in both unpolarized and polarized 
deep-inelastic scattering. The available resummed anomalous dimensions are discussed for all these cases, and 
the most important analytic and numerical results are compiled. The quantitative effects of these small- a; resum- 
mations on the evolution of the various parton densities and structure functions are presented, and their present 
uncertainties are investigated. An application to QED radiative corrections is given. 



1. Introduction 

The evolution kernels of both non-singlet and sin- 
glet, unpolarized and polarized parton densities 
contain large logarithmic contributions for small 
fractional momenta x. For unpolarized deep- 
inelastic scattering (DIS) processes the leading 
small- a; contributions in the singlet case behave 
like § 

where N is the Mellin variable. The leading terms 
for the unpolarized and polarized non-singlet and 
the polarized singlet cases are of the form |,| 

AT f a s\ k k 1 2fc-2 

The resummation of these terms to all orders in 
the strong coupling constant a s can be completely 
derived by means of perturbative QCD. Since in- 
finities, such as the ultraviolet and collinear diver- 
gencies, emerging in the calculation of the higher- 
order corrections have to be dealt with, the only 
appropriate framework for carrying out these re- 
summations is provided by the renormalization 
group equations. The impact of the resulting all- 
order anomalous dimensions on the behaviour of 
the DIS structure functions at small x depends as 
well on the non-perturbative input parton densi- 
ties at an initial scale Qq. Thus the resummation 
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effects can only be studied via the evolution over 
some range in Q 2 . 

This evolution moreover probes the anomalous 
dimensions also at medium and large values of x 
by the Mellin convolution with the parton densi- 
ties. Hence the small- a; dominance of the leading 
terms over contributions less singular as x — > in 
the anomalous dimensions does not necessarily 
imply the same situation for observable quanti- 
ties, such as the structure functions. These as- 
pects need to be considered to arrive at sound 
conclusions about the consequences of the small- 
x resummations on physical quantities. 

In the present paper we give a survey of the 
recent developments in this field. The general 
framework for the evolution of parton densities 
and structure functions is recalled in Section 2. 
Section 3 reviews the available results on the re- 
summed anomalous dimensions for the various 
DIS processes. Numerical coefficients for their 
expansions in a s are compiled, as well as the ana- 
lytical predictions for the most singular contribu- 
tions to the 3-loop splitting functions. The issue 
of subleading terms is discussed, guided by the 
known 2-loop results. In Section 4 the numeri- 
cal implication of these resummations are inves- 
tigated for the various unpolarized and polarized 
cases. The uncertainties due to possible less sin- 
gular terms and due to insufficiently constrained 
initial parton densities are illustrated. An appli- 
cation to QED radiative corrections is presented. 
Section 5 summarizes the main results. 
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2. The evolution equations 

The twist-2 contributions to any deep-inelastic 
scattering structure function can be represented 
in general, see e.g. ||,||, by the three flavour non- 
singlet combinations of quark densities 



9ns, 



Ins 



and the singlet quark and gluon distributions 



Qs 



N f 

i: = ^2(q r + q r ) 

r=l 



(1) 

(2) 
(3) 



Here Nf denotes the number of active (massless) 
quark flavours. The structure functions Fi(x, Q 2 ) 
are obtained by 



Fi{x,Q 2 ) 



2 Nf 

^ a ir Ci >r (x, Q 2 ) ® q r (x, Q 2 ) 



are identical for the combinations q^ s . and gj^jg. 
Since this is the case at all orders known presently 
(i.e. up to next-to- leading order, NLO) 2 , we will 
not investigate the evolution of the combination 
(||) separately in the following. 

The evolution equations for the non-singlet 
and singlet combinations of the parton distribu- 
tions are then given by 



dqj s (x,Q 2 ) 

d\nQ 2 
dq s {x,Q 2 ) 
d\nQ 2 



P^ s (x,a s )<g>q^ s (x,Q 2 ) , 
P s (x,a s )®q s (x,Q 2 ) . (8) 



The splitting functions P NS and Ps are speci- 
fied below. Note that the unpolarized (polarized) 
quark density combinations q^ s and g^ s evolve 
with P^ s (Pjjg) and P^ s (P^ s ), respectively. 

In the following, we will simplify the notation 
by dropping the subscripts 'NS' and 'S', and use 
the abbreviation a s = a s (Q 2 ) / '4ir for the running 
QCD coupling for convenience. The scale depen- 
dence of a s is governed by 



+ a lg c hg (x,Q 2 )® g{x,Q 2 ) , (4) 

where the factors depend on the electroweak 
couplings, and Cij(x,Q ) denote the respective 
coefficient functions. Finally ® stands for the 
Mcllin convolution in the first variable, 

A(x)®B(x) = (5) 
,i ,i 

dx\ I dx2 S(x — X1X2) A(x\)B{x2) ■ 
'0 Jo 

The above notation is used in a generic way for 
both unpolarized and polarized DIS, i.e. for the 
polarized case the replacements 



Ag , q^ Aq , and Ag, 



(6) 



are understood with, e.g., Ag given in terms of 
the spin projections gf and gj. via 



Ag = g? -gj. 



(7) 



Correspondingly the splitting functions (see be- 
low) and the coefficient functions in eq. ([|) have 
to be replaced. 

As long as the splitting functions P^ q and P^ 
(cf. ref. |Q) do not differ, the evolution equations 



da, 



tflnQ 2 



E« fc+2 

fc=0 



0k, 



(9) 



where only O = (11/3) C A - (A/3)T F N f and 
/3i = (34/3) C\ - (20/3) C A T F N f - AC F T F N f 
enter up to NLO. The colour factors are C F = 
(N 2 - 1)/(2JV C ) = 4/3, C A = N c = 3, T F = 1/2. 
For the numerical calculations in Section 4 we use 



a s = -^ln(Q 2 /A 2 ) 



1 - 



/?! lnln(Q 2 /A 2 ) 



P 2 HQ 2 /A 2 



(10) 



with A being the QCD scale parameter. The 
splitting functions and coefficient functions can 
be represented by the series 



P ± (x,a s ) 
P(x,a s ) 



Pr~(x) , 



(11) 



2 So far only for the anomalous dimension 7^0 the first 
moments have been calculated to 3-loop order M . 
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i(x,Q 2 ) = 6(l-x)6 jq + 



^2a l s c ijtl (x) 
1=1 



(12) 



Unless another scheme is stated explicitly, we will 
always refer to the MS scheme both for renor- 
malization and factorization, and take Q 2 as the 
renormalization and factorization scale. 

The expansion coefficients Pf{x) and P" npol (x) 
are subject to the sum rules 



dxPf(x) = 0. 



dx . 



(13) 



which are due to fermion number and energy mo- 
mentum conservation, respectively. By now all 
the unpolarized and polarized splitting functions 
are completely known up to NLO, 1 = 1. The full 
expressions for their x-dependences can be found 
in refs. fl— The most singular c ontr ibutions 
as x — > will be displayed in Section \5A\ 

The parton densities are not observables be- 
yond leading order. Hence it is convenient 
to consider also the evolution equations for re- 
lated physical quantities, the structure functions 
Fi(x,Q 2 ), directly. In non-singlet cases, the all- 
order resummation of the leading small- a; contri- 
butions has in fact been given 0] on the level of 



the structure function combinations F; 
Their evolution equations read 

1 



: (x,Q 2 ). 



dF±(x,a s 



da s f3 a 2 

after a transformation to an equation in a 
e.g., the NLO kernels can be written as 

KfAx,a s ) = 



F*(x,a s )(U) 
Here. 



(15) 



a s P (x) +a 2 s 



Pt{x)-^P {x)-l3 Q ct <1 {x) 



with c i (x) denoting the corresponding coeffi- 
cient function combinations. Generally the terms 
cx a s (a s In 2 x) 1 emerge in the a s expansion of 
the kernels K^(x,a s ) only in combination with 
the coefficient (3q. As will be outlined below the 
leading small- a: contributions to these kernels co- 
incide with those of the MS splitting functions 
Pj^ix), since the corresponding coefficient func- 



tions Cj ; (x) turn out to be less singular as x 



0. 



3. Resummation of leading small- a; terms 

3.1. The non-singlet case 

The most singular contributions to the Mellin 
transforms of the all-order evolution kernels 
K {x, a s ) for the non-singlet structure functions 
have been obtained via 

M [K^ (a s )] (N) ee f 1 dx x N - 1 K^ Q (x 1 a s ) 
Jo 



8tt 2 



f£(N,a 3 ) (16) 



from positive and negative signature amplitudes 
f^(N,a s ) studied in ref. @ (cf. Section 3.2): 




Here the colour octet amplitude f^(N,a s ) reads 

f+(N,a s ) = 16ir 2 N c a s ^ln [e z2 ^D_ 1/[2NS] (z) 

(18) 



with z = N/\/2N c a s and Dy(z) denoting the 
parabolic cylinder function |12| . 

Expanding the resummed anomalous dimen- 
sions (17) into a series in a s and transforming 



I 




K i 





2.667E0 


2.667E0 


1 


3.556E0 


5.333E0 


2 


1.580E0 


1.432E0 


3 


3.512E-1 


9.964E-1 


4 


4.682E-2 


-2.078E-1 


5 


4.162E-3 


1.448E-1 


6 


2.643E-4 


-5.777E-2 


7 


1.258E-5 


2.168E-2 


8 


4.661E-7 


-7.173E-3 


9 


1.381E-8 


2.143E-3 


10 


3.348E-10 


-5.827E-4 



Table 1: The coefficients K, of the expansion of 
K^ Q (x,a s ) in terms of a s (a s In x) as obtained 
from the resummations in eqs. ([[7]) and (|j~8|). 
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back to a;-space yields the numerical coefficients 
shown in Table 1. The first two terms of the re- 
summed kernel K^_^ (x,a s ) agree with the lead- 
ing small- x contributions of the corresponding LO 
and NLO splitting functions P ± (x,a s ) @. This 
expansion also reveals a significant theoretical dif- 
ference between T^ Q (N,a s ) and r~_^ (iV, a s ): 
the a s series is convergent in the former but not 
in the latter case, where it involves the asymp- 
totic expansion of D p {z). 

Unlike the splitting functions, the coefficient 
functions cf(x, a s ) are known up to next-to-next- 
to-leading order (NNLO), I = 2, in the MS 
scheme |l4|,[l5|]. At small x they rise only as 

(19) 



(',"-, a Ins , cf 2 oc hi* x 



Therefore the third expansion coefficient of 
K^_, (x,a s ) leads to a prediction for the most 
singular parts of the non-singlet MS splitting 
functions in NNLO, P 2 ± (x), given by [pf 



^2, x-*o{ x ) 



10 



C% + 4C f Ca — CfC A 



(20) 



ln 4 x. 



All the methods of this section have been applied 
to QED and also comparisons with the available 
fixed order calculations were carried out. For de- 
tails we refer to refs. p||l7| |. 

3.2. The polarized singlet case 

The amplitude relations of ref. [g] have been gen- 
eralized to the polarized singlet case recently Pj : 



F {N,a s ) = 16tt 2 ^M 



N 



(21) 



N 2 



F 8 (iV,a s ) = 167r 2 ^M 8 
2a s „ d 



(22) 



1 1 



-— C G — F 8 (N,a s ) + -^r—F 2 JN, a s ) 
N G dN 8V ' sj 8tt 2 N 8V ' ' 



Here the basic matrices are given by 



Mo - | 2Cp Wa I , Go - 



(C F -C A /2 —TpNf 



C F 

o c A 



C A 



2C A 



(23) 



Note that the matrix Mq is the x — > limit of 
the well-known matrix of the polarized splitting 
functions in LO B . The equations for (N, a s ) 
and fg(N, a s ) of Section 3.1 are entailed in these 
expressions by keeping only the qq-entvies of the 
matrices (p3|), and, in the '+'-case, additionally 
dropping the -Ps-term in eq. (21). Also for the 
polarized singlet structure function gi(x, Q 2 ) the 
coefficient functions are known up to NNLO Q, 
and their leading small- a; behaviour is the same 
as in eq. (^]). The resummed leading small- ie 
contributions to the splitting functions are thus 
related to F a (N,a s ) by 



(24) 



2=0 



—M-^FoiN,^)] (x) 



Eqs. ( ^l| ) and (^3|) have been solved directly in 
terms of a series in a s in ref. |l8fl . Unlike the pre- 
vious non-singlet case, the representation is 
needed for the analytical iV-space solution of the 
evolution equations here, cf. Section 4.3. Also in 
the present case the lowest-order expansion co- 
efficients PSa°{x) agree with the corresponding 
limit of the NLO splitting function matrices 0. 
The predictions for the NNLO quantities P^^? 
read 111: 



R 



x — *0 
99,2 



5C 



, - - C A + 6 Ca Cf 



In x 



p x^O 



(x) = 



- 8 TpNf Cf — 6T F Nf C A 

15C A + 2C 2 F -6C F C a 



\t f N s 



\T F N f C F 



In 4 x 



(25) 



Pg q a\x) =-C F 15C 2 a-2C 2 f + 6CfCa 



8T F N f C F 



In x 



p x^O 



(x) 



2 r 
3 



28 C\ + 2 T F N f C 2 A - AT F N f C 2 F 
24C F TpNf Ca] ln 4 x. 



The expansion coefficients up to I = 10 are shown 
in a compact numerical form in Table 2. For a 



more precise representation sec ref. 18 1. 
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1 
2 
3 
4 
5 
6 
7 
8 
9 
10 





1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



N f =3 



AD 



2.667E0 
-1.067E1 
-3.679E1 
-6.642E1 
-6.110E1 
-3.858E1 
-1.680E1 
-5.632E0 
-1.424E0 
-2.991E-1 
-4.869E-1 



AD 



-6.000E0 
-4.400E1 
-1.394E2 
-2.288E2 
-2.154E2 
-1.347E2 
-5.955E1 
-1.979E1 
-5.082E0 
-1.050EO 
-1.757E-1 



,(0 



5.333E0 

3.911E1 

1.240E2 

2.034E2 

1.915E2 

1.197E2 

5.294E1 

1.759E1 

4.517E0 

9.331E-1 

1.562E-1 



AD 



2.400E1 
1.280E2 
4.189E2 
6.981E2 
6.685E2 
4.201E2 
1.868E2 
6.213E1 
1.602E1 
3.303E0 
5.557E-1 



N f =4 



A') 



2.667E0 
-1.600E1 
-4.953E1 
-8.573E1 
-7.633E1 
-4.649E1 
-1.956E1 
-6.326E0 
-1.546E0 
-3.133E-1 
-4.923E-2 



An 



-8.000E0 
-5.867E1 
-1.788E2 
-2.859E2 
-2.595E2 
-1.569E2 
-6.688E1 
-2.148E1 
-5.319E0 
-1.063E0 
-1.713E-1 



,(0 



5.333E0 

3.911E1 

1.192E2 

1.906E2 

1.730E2 

1.046E2 

4.459E1 

1.432E1 

3.546E0 

7.086E-1 

1.142E-1 



An 



2.400E1 
1.227E2 
3.905E2 
6.316E2 
5.831E2 
3.540E2 
1.519E2 
4.882E1 
1.215E1 
2.421E0 
3.928E-1 



Table 2: The coefficients P {[) of the expansion 
of P(x, a s ) x ^o in eq. ( [24] ) in terms of a s (a s In 2 x) 1 
for three and four quark flavours. 

The leading small- a; off-diagonal elements of 
P(x, a s ) x —>o are related by 

P q x g 7°(x)/(T F N f ) 



(26) 



for all I. Also the case of an Af — 1 supersym- 
metric Yang-Mills field theory i.e. Ca — Cf = 1, 
Nf — 1, and Tp = 1/2, has been considered. The 
so-called supersymmetric relation 

P qq ,l(x) + PmM ~ PlM ~ PggM = (27) 

is satisfied for the small-a: leading terms. In fact 
even more restrictive relations are fulfilled in this 
case, cf. ref. |l|]. Finally it should be noted that 
there is no overlap of the present small- a; resum- 
mation with the large-iVf expansion Jl9| of the 
all-order splitting function matrix. 



3.3. The unpolarized singlet case 

Unlike the cases discussed in the previous sec- 
tions, where the leading small- a; singularity in the 
complex N plane is situated at N = 0, the corre- 
sponding poles of the anomalous dimensions for 
the unpolarized singlet evolution are located at 
N = 1. The all-order resummation of the most 
singular contributions as x ~ * 0, 7l(^V, a s ), to 
the anomalous dimensions was derived in It 
is found as the solution of 



N = AC A a sX h L {N,a s )] , 
with 

X ( 7 ) = 2^(1) - V(7) - V>(1 - 7) 



(28) 



(29) 



Here ip(z) denotes the logarithmic derivative of 
Euler's T-function. 7^ is a multi-valued function 
for complex N. The perturbative branch is se- 
lected by requiring 



"f L {N,a s ) 



4:C A a s 
N-l 



for \N\ 



(30) 



when solving eq. (|2§|). The singularity struc- 
ture of the solution in the complex N plane was 
studied in detail in refs. |2^,|2l|. Finally the re- 
summed small- a; contributions ~/l(N) to the sin- 
glet anomalous dimension matrix ~f(N), related 
to the splitting functions in eq. (11) by 



7(JV, o.) = -2 f dx x N ~ 1 P{x, a.) , (31) 
Jo 

are in this approximation obtained by 



7lW as) 





Cf/Ca 



j L (N,a s ). (32) 



The 0(a s ) correction 7 NL (iV, a s ) to this most 
singular part as x — > have been calculated 
in ref. Q for the quark anomalous dimensions 
"fqq,qg- ^ n ^ ne ^IS factorization scheme the ma- 
trix of these next-to-leading small- a; contributions 
is given by 



7 NL (iV,a s ) = -2 



'NL — 3C 
7gg,NL 



sTf] 7NL 
7s9,NL 



(33) 
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where 7 nl is used as an abbreviation for the func- 
tion 7^(iV, a s ). It can be recursively expressed 
by 7i(JV, a 3 ) and reads 



2AT F a. 



2 + 3 7L -3 7 £ [B{l- lL ,l + lL )f 



3-2 7L P(2 + 2 7l ,2-2 7l , 
with B(x,y) denoting the Beta function and 

r(i- 7 )x(7) 11/2 



(34) 



r i 



exp 



7){- 
7^(1) + 



7X'(7)} 
o 



(35) 



V-'(i-C) 



x(C) 



The calculation of the yet unknown gluonic en- 
tries 7 sg ,NL and 7 99 ,nl in eq. (||) is in progress 
]22] , p3] |. A first contribution to 7 9g ,NL oc Nf 
has been determined recently [^3| in the Qq- 
scheme |^4|,^5| , which has been introduced in the 
framework of /c^-factorization. 

Both 7l and 7 NL can be represented by in- 
finite series in a s . The analytic expressions are 
straightforwardly obtained but are rather lengthy. 
The numerical size of the resulting coefficients A\ 
and Bi in the DIS scheme is illustrated in Table 3. 
The matrix P(x,a s ) Tes of the splitting functions 
including the small-! resummed terms (cf. Sec- 
tion 4.4) beyond the complete LO and NLO ma- 
trices from the fixed-order calculations reads: 



P(x, a s )°f = a s P (x) + aiP x {x) 



DIS 



1=2 



■ In' 





C F IC A 1 



(36) 



i=i 



:5io +a -ln' - 



Cf/Ca 




In the subsequent numerical treatment we will use 
the labels 'Lx ' for results obtained with the lead- 
ing series {A{), and 'NLx' in the cases where the 
Pi-terms in eq. (36) have been taken into account 
additionally. 

Finally we list the presently available predic- 
tions from the small- x resummations for the most 
singular contributions P^J of the NNLO split- 
ting functions. There are no afln 2 x terms, due 
to the matrix structure of eq. (|32|) and the van- 
ishing of A2 in eq. (p6[). For the quark splitting 



I 


A l 


Bi 





1.200E1 


N r 


3.467E1 


1 


0.000E0 


Nf 


4.415E2 


2 


0.000E0 


Nf 


9.528E3 


3 


8.308E3 


Nf 


6.877E4 


4 


0.000E0 


Nf 


4.040E5 


5 


5.160E4 


Nf 


3.411E6 


6 


8.629E5 


Nf 


1.293E7 


7 


1.721E5 


Nf 


5.518E7 


8 


5.104E6 


Nf 


2.601E8 


9 


2.879E7 


Nf 


7.086E8 


10 


1.433E7 


Nf 


2.343E9 



Table 3: The coefficients Ai and Bi of the expan- 
sion for the small-a; resummed splitting functions 
in the DIS scheme, see eq. (^6|). For completeness 
also Aq, Ai and Bq are given. For a more precise 
representation and higher-Z terms, cf. refs. |2g , |27]| . 

functions || the a 3 s In x terms read in the DIS 
scheme 



px-fOf„\T>IS 



«?,2 



568 



8 



C F C A -\n( - 

x \x 



The corresponding MS results are given by 

r 



9 x 



??,2 

and 



(37) 



(38) 



(39) 



in both schemes. Unlike the cases discussed in the 
previous sections, the coefficient functions con- 
tain terms as singular as the splitting functions 
in the MS scheme. 

3.4. Less singular contributions 

The terms in the splitting functions Py and P ± , 
which are less singular by one (or more) powers of 
ln(l/x) as x — > than the leading contributions 
discussed in the previous sections, are presently 
unknown in almost all cases. Such subleading 
contributions, however, can potentially prove to 
be as important as the leading terms, since the 
splitting functions and coefficient functions enter 
observable quantities always via Mellin convolu- 
tions with the parton distributions. 
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This situation can be illustrated by a simple 
example, cf. ref. ^4j. Consider the lowest-order 
gluonic contribution to the longitudinal structure 
function Fl(x, Q 2 ), given by 



(40) 



[ ld y y l(l-y)*g( X /y,Q>) 

Jx v y 



If one replaces the term 1 — y originating from the 
coefficient function by its small-?/ approximation 
1 for small values of x, the result for F£ changes 
by a factor of about 4 for typical parametriza- 
tions of the gluon density! Due to the Mcllin 
convolution and the fact that g{x) becomes very 
large as x — ► 0, the coefheient function at medium 
and large y contributes essentially. On the other 
hand the coefficient function in the range y^x, 
where the small- y approximation is justified, sam- 
ples g(x/y<l) which is however small. Similar 
observations can be made for other convolutions 
considered in the previous sections as well. 

The non-singlet '-' and the unpolarized singlet 
splitting functions are constrained by conserva- 
tion laws, see eq. (|l3|). The resummed contribu- 
tions discussed in Sections 3.1 - 3.3 do not obey 
these constraints, however, less singular terms re- 
store these sum rules. Also for the cases in which 
the anomalous dimensions are not subject to such 
constraints, less singular terms with sizeable co- 
efficients exist for example in NLO, e.g. eq. (pLl]). 

In order to evaluate the possible impact of such 
terms, their numerical coefficients have to be esti- 
mated. At present the only source of information 
are the fully known LO and NLO splitting func- 
tions. The dominant and subdominant terms as 
x — > for the NLO anomalous dimensions are re- 
called in eqs. (f4l|)-([43"|). The results are presented 
in the general form, as well as, for easier compar- 
ison of the numerical size of the coefficients, in- 
serting the number of active flavours as used in 
the numerical applications in Section 4. In the 
non-singlet cases one finds in the MS scheme 
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Note that the subleading terms are roughly of the 
same size in both cases, despite only one of the 
combinations being constrained by a sum rule. 
The corresponding results for the polarized sin- 
glet case, also in the MS scheme, are given by 
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The first two terms of the unpolarized singlet 
anomalous dimensions expanded at N = 1 read 
in the DIS scheme 

t£S = ^^(y + 2.56 + 100.83/V/ 

7Vf=4 123.3 
~ N- 1 

7 9 Ti W-o = 3=|^^ + 218.67JV} - 1.787V/ 2 
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~ N - 1 



846.2, 
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One notices that the subleading terms occur in 
general with signs opposite to those of the domi- 
nant ones. Their prefactors are of the same order, 
but in most cases a factor of about 2 to 4 larger. 

Thus introducing subleading terms with pref- 
actors up to two times larger than those of 
the leading terms appears to yield reasonable 
and conservative estimates for the possible im- 
pact of subleading terms. The following modifi- 
cations of the resummcd anomalous dimensions 
T(N, a s ) have accordingly been studied within 
refs. ||l|BM8p|B7l: 



A 
B 
C 
D 



T(N, a.) 
T(N, a,) 
T(N, a.) 
T(N, a.) 



T(N, a s )-T(l,a s ) 
T(N, a s )(l-N) 
T(N, a s )(l - 2N + N 2 ) 
T(N, a s )(l-2N + N 3 ) 



(44) 



where the replacement N — ► N — 1 is understood 
for the case of Section 3.3. 

Let us finally discuss also the case of the '-' 
non-singlet evolution in QED. For the evolution 
kernel also the terms of 0(a 2 \nx) were calcu- 
lated for e + e~ annihilation in the on-mass-shell 
scheme (OMS) in ref. pi]. All contributions but 
those due to the vacuum polarization diagram 
cancel in this order. One obtains 
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In 2 x 
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N 3 
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(45) 



Unlike for most of the examples discussed above, 
in this particular case the term being suppressed 
by one order in In x has thus a smaller coefficient 
than the leading singular contribution. 

4. Numerical consequences 

4.1. The unpolarized non-singlet case 

The evolution equations (Q) for the non-singlet 
combinations of structure functions can be solved 



analytically in Mellin-A space. Taking into ac- 
count the resummed kernels Kjjr_^ Q of eq. (|l^) in 
addition to the full splitting functions up to NLO, 
the solution can be written as [Jl3.l6| 
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do 
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(46) 



with 



dxx N - 1 Pi(x) 



dxx Ci,l(x) , 



(47) 



and ao = a s (<3o)- Here ~ff{N) denotes the two- 
loop anomalous dimension ^(N) with the lead- 
ing I/TV 3 term subtracted. The effect of this term 
is included to all orders in the exponential factor, 
which in turn is obtained from eq. ( |l7j ) by remov- 
ing the LO contribution included in 70 (N): 



r ± (Af, a s ) 



T^ (N,a s ) 
T^ (N,a s ) 



— lim [Nj (N)} 



4C F 



(48) 



N 

The NLO evolution of F~-(N, a s ) can be recov- 
ered from eq. ( |46"1) by expanding the exponential 
to first order in a s and ao. The inverse Mellin 
transformation of the final results back to x-space 
is performed by a numerical integral in the com- 
plex iV-plane, see e.g. ref. pj| . 

The remaining quadrature in (|4^) can be per- 
formed analytically for the '+'-case pq] , and has 
to be done numerically for the '-'-combinations 
involving the parabolic cylinder function D p (z — 
N/y/2N c a s ). Additional information can be 
obtained by expanding the resummed kernels 
r + (A^, a s ) and r~(A^, a s ) in a s , in the latter case 
using the asymptotic expansion of D p (z) p2] |. 

The evolution of the '-'-combination 



xF»(x, Ql) = \ lxF£ N {x, Ql) + xF* N 



2 L"3 

= Cp (x,Qq)<8 [xu v + xd v ] (x, 



0) (49) 



9 



for an isoscalar target N and the '+ '-combination 



F 2 ep 



pen 



(x,> 



} (x, Qq) ® - [xu v — xd v — 2(xd — xu)] (x, 



(50) 



have been investigated in refs. ]l^,[l6|. As in all 
other numerical examples displayed below, the 
reference scale for the evolution (46) is chosen as 
Qq = 4 GeV 2 , and the same input parameters are 
employed for the NLO and the resummed calcu- 
lations. In the present case, the initial parton dis- 
tributions have been adopted from the MRS (A) 
global fit |}0| together with the value of the QCD 
scale parameter, A^(Nf — 4) = 230 MeV. 
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Figure 2: The same as in Figure 1, but for the 



structure function combination F 2 ep - 



F 2 en . Instead 



of the prescription '(A)', the result without any sub- 
leading terms is shown for this '+'-case. 

are rather 'steep': their rightmost singularities in 
the complex iV-plane lie about 0.5 units or more 
to the right of the leading singularity of the non- 
singlet anomalous dimensions at N = 0. 

In Figures 1 and 2 the NLO results for xF^ 
and F 2 P — F 2 n are shown, together with the corre- 
sponding resummation corrections, down to x as 
low as x = 10~ 15 . Even at these extremely small 
values of x, the effect of the resummed anoma- 
lous dimensions stays at the level of 1% or be- 
low. This is in striking contrast to the expecta- 
tion of ref . J31J , where corrections of up to factors 
of 10 in the HERA kinematical regime were an- 
ticipated. Moreover the resummation effects re- 
main very sensitive to presently unknown terms 
less singular as x — > in the splitting functions. 
This is illustrated by the impact of the prescrip- 
tion (B) of eq. (0), which removes as much as 
about two thirds of the resummation effect even 
at asymptotically low x. 

Another interesting issue is the a s -expansion of 
the resummed kernel ( |l7| ) mentioned be- 

xd v (x, Qq) ~ x 33 1 30 1 . Hence these distributions fore. In Figure 3 the resummation corrections 



Figure 1: The small-x Q 2 -evolution of the non- 
singlet isoscalar structure function xF 3 N in NLO and 
the absolute corrections to these results due to the 
resummed kernels of Section 3.1. '(A)' and '(B)' 
denote two prescriptions for implementing fermion 
number conservation, see eq. (p4|). 



The small- a; behaviour of the most relevant initial 

and 



distributions is given by xu v (x,Qq) 



„0.54 
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with T ± (N,a s )/a s in eq. ( fl6| ) expanded to order 
I are compared to the full results for the cases 
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Figure 3: The ratio of the resummation cor- 
rections with the resummed anomalous dimensions 



- x— >0 



expanded at order a l s +1 to the complete ef- 
/ F 3 N (upper part 
F 2 en (lower part). 



fects for the '-'-quantity F 3 N (upper part) and the 



-combination F 2 ep 



The asymptotic series for the '-'-case leads to a 
good approximation at x > 10 -6 , but starts to di- 
verge below. In the '+'-case, on the other hand, 
the Taylor series converges in the whole x-range 
considered. In both cases the next two terms be- 
yond NLO, I = 3, contribute more than 90% of 
the final resummation effect, again even down to 
x = 10~ 15 . The reduced stability for F 2 ep - F 2 en 
around x = 10~ 7 is immaterial, since A/ changes 
sign here. 



4.2. The polarized non singlet case 

The solution of the evolution equations proceeds 
in the same way as in the previous section, see 
eqs. (fl6|)-(48). The present case is however prac- 
tically even more interesting. Firstly the non- 
singlet structure functions are, unlike in the unpo- 
larized case, not a priori suppressed versus their 
singlet counterparts at very low x. Secondly the 
shapes of the polarized parton densities are not 
well established yet ||3| by the experimental re- 
sults. This is illustrated in Figure 4, where the 
small- x extrapolations of the structure function 
gl P — g{ n (see eq. (^l|) below) are compared for 
two choices of the initial distributions. 
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Figure 4: The NLO evolution of the polarized 
non-singlet structure function combination g± p — 
gf n for the input densities from refs. [ fP] ] and 



At small x the valence quark densities of CW 
Q are relatively flat: Au v , Ad v (x, 10 GeV 2 ) ~ 
2,-0.17, +0.29^ rp ne more recen ^ distributions of 

GRSV |55| are on the other hand approxi- 
mately as steep as the unpolarized distributions, 

,-0.67 



and Ad v (x, Qq) 



For the evolution of the former (latter) input 
A m (N f = 4) = 230 (200) MeV is employed in 
this section. 



11 



The effect of the resummed anomalous di- 
mensions (|l7]) has been investigated in 
refs. M,[l6i for the '-'-combination 



10 



9i P {x,Ql)-gr{x,Ql) 



(51) 



Cg ± (x, Ql) ® - [Au„-Ad„ + 2(Au-Ad)] (x, Q%) . 

Here also the '+'-case of the 7Z- interference 
structure function is considered: 

a^yzfa Qo) = 4A X > Ql) ® AdJ(x, Qo) • 

(52) 

The resummation corrections for g* p — gf n and 
9b P -yZ are sri0wn m Figures 5 and 6, respectively. 
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Figure 5: Relative corrections to the NLO small-x 
Q 2 -evolution of the '-'-combination g^ p — gf n due 
to the resummed kernel of Section 3.1 for the initial 
distributions of refs. Q and ||5). 'A', 'B', and 
'D' denote different prescriptions for implementing 
fermion number conservation, see eq. (|44 
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Figure 6: The small-.! evolution of the non-singlet 
interference structure function g^ P lZ in NLO and 
the resummation corrections to these results. The 
possible importance of less singular terms in the 
higher-order splitting functions is illustrated by the 
prescription '(B)' for this ' + '-case. 

For the relatively flat CW input p4| ], the re- 
summation effect on g^ v — gf n reaches 15% at 
x = 10~ J . However, in the restricted kinemat- 
ical range accessible in possible future polarized 
electron-polarized proton collider experiments at 
HERA fl37|] , it amounts to only 1% or less. For the 
steeper GRSV initial distributions , the effect 
is of order 1% or smaller in the whole x range, as 
in the unpolarized cases considered above. Hence 
also here the results do not come up to previous 
expectations of huge corrections up to factors of 
10 or larger jj2| in the HERA range. The out- 
come is very similar for g§ P lZ , see Figure 6. 

As for the unpolarized non-singlet structure 
functions, the resummation results also in the po- 
larized cases are not stable against possible sub- 
leading contributions in the higher-order anoma- 
lous dimensions. With respect to the a s expan- 
sions of the kernels (17) the situation is 



also similar to the unpolarized cases. 
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4.3. The polarized singlet case 

The solution for the singlet part of the evolution 
equation (||) cannot be given in a closed form be- 
yond leading order. This difference to the non- 
singlet case is due to the non-commutativity of 
the splitting functions matrices Pi (x) in eq. ( pd| ) 
for different orders in the strong coupling a s . 
Thus the solution has to be written down as a 
power series in a s , in iV-space resulting in 

q (N,a s ) = [l + E-^Hfe) 70 ™ 

l + Y,4Ui(N)Y 1 q(N,a ). (53) 
i=i 

Here ~f (N) is related to the matrix of the LO 
splitting function Pq(x) as in eq. (|47|), and as be- 
fore ao = o, s (Qq). The singlet evolution matrices 
Ui(N) can be expressed in terms of the anoma- 
lous dimensions Jk<i(N)- Technical details can 
be found in rcf. [ p7[ . 
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Figure 7: The evolution of the polarized singlet 
combination xAY, in NLO and including the re- 
summed kernels. The impact of possible sublead- 
ing terms is illustrated by the prescription '(B)' in 
eq. (0). The input densities are from ref. 



The numerical consequences of the resumma- 
tion in Section 3.2 on the polarized parton densi- 
ties and the structure functions g^ p and gf n have 
been investigated in ref. p8[. In Figures 7 and 
8 the results are displayed for the polarized sin- 
glet and gluon densities, xAYi and xAg, respec- 
tively. The initial distributions at Qq = 4 GeV 2 
have been adopted from the GRSV 'standard' 
parametrization |55| . The evolution has been per- 
formed for A m (N f = 4) = 200 MeV, and - dif- 
ferent from all other cases shown in this paper - 
with only three active quark flavours in the split- 
, i.e. AE is given by 



ting functions |35 3 
AE = Au 



Au + Ad + Ad + As + At 



(54) 



As in the corresponding non-singlet '-'-case, the 
expanded solution ([53]) represents an asymptotic 
series, which diverges at very small values of 
x. For all the experimentally relevant cases, 
x > 10~ 5 , retaining 8-10 terms in eq. ([33]) is how- 
ever adequate for obtaining accurate results. 
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Figure 8: As in Figure 7, but for the polarized 
gluon momentum distribution xAg. As in the pre- 
vious figure, the Q 2 -values in the legend are ordered 
according to the sequence of the curves at small x. 
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Figures 7 and 8 show that the resummation 
effects are much iarger for AS and Ag than for 
the non-singlet quantities considered in Section 
4.2, as to be expected from the comparison of the 
expansion coefficients in Tables 1 and 2. E.g., the 
ratio of the (unsubtracted) resummed results to 
the NLO evolution amounts to about 1.72 (1.64) 
for AS (Ag) at Q 2 = 10 GeV 2 and x = 10~ 4 . 

Also illustrated in these figures [by the results 
for the prescription '(B)' in eq. (|44|)] is the pos- 
sible impact of the yet uncalculated terms in 
the higher-order anomalous dimensions which are 
down by one power of N with respect to the re- 
summed leading pieces as N — > 0. The effect of 
these additional terms can be very large, in the 
present example the resummation correction be- 
yond NLO is practically cancelled. 



Q 2 


10 GeV 2 


100 GeV 2 


X 


10- 4 


10~ 3 


io~ 4 


10" 3 


xAS 


-0.0100 
-0.0285 
-0.0473 


-0.0169 
-0.0396 
-0.0560 


-0.0171 
-0.0505 
-0.0855 


-0.0218 
-0.0523 
-0.0772 


xAg 


0.019 
0.101 
0.201 


0.034 
0.152 
0.294 


0.053 
0.226 
0.432 


0.071 
0.281 
0.528 



Table 3: A comparison of the resummed evolu- 
tion of the polarized parton distributions for dif- 
ferent assumptions on the gluon distribution Ag. 
Upper lines: minimal gluon, middle lines: standard 
set, lower lines: maximal gluon (and corresponding 
quark distributions) of ref. J||] at Q 2 Q = 4 GeV 2 . 

The small- a; evolution depends strongly also 
on the virtually unknown [|33"| gluon input den- 
sity. This is obvious from Table 3, where the re- 
summed results of Figures 7 and 8 are compared 
at two representative values of x and Q 2 to those 
obtained by evolving in the same way the 'mini- 
mal A<?' and 'maximal Ag' distributions of GRSV 
|35| . The variations are up to a factor of almost 
5 (10) for AS (Ag), respectively. Thus both the 
unknown less singular terms in the anomalous di- 
mensions and the present bounds on Ag, which 
are rather weak still, are the dominant sources of 
uncertainty at small x. 



4.4. The unpolarized singlet case 

The solution of the evolution equations for the 
unpolarized singlet parton densities is analogous 
to the polarized case considered in the previous 
section, see eq. (|53|). The present case is special - 
and has therefore attracted much interest over the 
past years - since only here precise measurements 
have been performed for small x, at HERA pp| . 
The quantitative impact of the resummation dis- 
cussed in Section 3.3 has been studied for parton 
distributions and structure functions in refs. |2o| 
and [^7j . The latter analysis confirms and extends 
the former one. Related investigations have been 
carried out in refs. |2^ , |39| . 

Below the results are shown for initial distribu- 
tions which, although representing a somewhat 
simplified input, incorporate all features relevant 
to this study in a sufficiently realistic way, espe- 
cially the small- a; powers as supported by HERA 
structure function data |^(| . Specifically, we take 
in the DIS factorization scheme at Q\ — 4 GeV 2 : 

xu v = A u x°- 5 (1 - x) 3 , xd v = A d x - 5 , (1 - x) 4 
xS = S — xu v — xd v — Asx~ ' 2 , (1 — x) 7 
xg = A g x~°' 2 (l — x) 5 , xc = xc = 0. 

(55) 

The evolution is performed for four active (mass- 
less) flavours, using A-^(N f = 4) = 250 MeV. 
The (SU(3)-symmetric) sea is assumed to carry 
15% of the proton's momentum at the input scale; 
together with the sum rules this fixes the prefac- 
tors Ai. 

Figures 9 and 10 compare the resummation re- 
sults separately for the Lx Jl]] and NLx || series 
(c.f. Section 3.3) to the standard NLO evolution. 
In the Lx case, as expected from the matrix struc- 
ture in eq. ([32]), the main effect is exerted on the 
gluon density xg. The impact on the quark evo- 
lution is rather moderate. Specifically, the ratio 
to the NLO results amounts to about 1.3 (1.03) 
for xg (xS), respectively, at Q 2 = 10 GeV 2 and 
x — 10~ 4 , taking the prescription '(A)' of eq. ( |44| ) 
for restoring the energy-momentum sum rule. In- 
cluding also the NLx quark terms || , on the other 
hand, results only in a small further modification 
of the gluon evolution, whereas the quark distri- 
butions are drastically affected. The ratios to the 
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Figure 9: The small-x evolution of the total unpo- 
larized sea quark density xS including the resummed 
Lx Jlj] and NLx kernels H as compared to the 
NLO results. Two prescriptions for implementing 
the energy-momentum sum rule have been applied, 
c.f. Section 3.4. 

NLO results now read 1.3 (3.1) for xg (xS) under 
the same conditions as before. 

A flavour of the possible importance of 
presently unknown less singular terms in the 
higher-order anomalous dimensions is provided 
by the difference of the results of the choices ' (A) ' 
and '(D)' in eq. (Q). Such terms can be vitally 
important, like in the polarized case studied in 
the previous section. As obvious from the figures 
not even the sign of the deviation from the NLO 
evolution can be taken for granted. 
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Figure 10: As in Figure 9, but for the unpolarized 
gluon momentum distribution xg. 

4.5. QED non singlet radiative corrections 

The resummation of the 0(a In 2 x) terms may 
also yield non-negligible contributions to QED 
radiative corrections jl6). This has been investi- 
gated recently for the case of initial-state radia- 
tion in deep inelastic eiV scattering. In the range 
of large y the effect can reach around 10% of the 
differential Born cross section |{[7| , see Figure 11. 
These terms are not covered by the higher order 
resummations studied so far ]4l| , ^2| and reduce 
their effect JlTj. Yet the complete NLO correc- 
tions for this process are not known and the size 
of less singular terms at 0(a 2 ) and their impact 
on the QED corrections is still to be determined 
(see, however, eq. (fl5|)). 
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Figure 11: The non-singlet (NS) resummed con- 
tribution of the 0(a(a In 2 x) 1 ) terms for the QED 
initial-state correction to the neutral current deep- 
inelastic scattering cross section at HERA. The re- 
sults are shown normalized to the differential Born 
cross section. 

For the corresponding corrections to cr(e + e _ — > 
near the Z-pe&k the situation is different. 
Since there the QED correction are completely 
known up to 0(a 2 ) ^8|, a correction due to the 
resummation of the 0{a l+1 ln 2i x) terms is only 
necessary for the terms beyond two loop order, as 
in the QCD cases considered before. A numerical 
study of these effects can also be found in ref . |l7j . 

5. Summary 

The resummations of the leading small- a; terms 
in both unpolarized and polarized, non-singlet 
and singlet anomalous dimensions have been dis- 
cussed. At NLO the results agree with those 
found for the most singular terms as x — > 
in fixed-order calculations. The so-called su- 
persymmetric relation is satisfied by the results 
for the most singular small- a; terms to all or- 
ders, again for both the unpolarized and polar- 
ized cases. These resummations allow the predic- 
tion of the leading x — > contributions to the 



3-loop (NNLO) anomalous dimensions in the MS 
scheme |||^,|l6|,|l8| . The coefficient functions are 
less singular for the non-singlet and polarized sin- 
glet cases up to NNLO, 0{a 2 s ). 

For the non-singlet structure functions the cor- 
rections due to the a s (a s In 2 x) 1 contributions are 
about 1% or smaller in the kinematical ranges 
probed so far and possibly accessible at HERA 
including polarization | p^| , ^6[ . The non-singlet 
QED corrections in deep-inelastic scattering re- 
summing the 0(a In 2 x) terms can reach values of 
up to 10% at x fa 1(T 4 and j>0.9@. 

In the singlet case very large corrections are ob- 
tained for both unpolarized and polarized parton 
densities and structure functions [|l8|,^o|,^7| . As in 
the non-singlet cases possible less singular terms 
in the higher order anomalous dimensions, how- 
ever, are hardly suppressed against the presently 
resummed leading terms in the evolution: even 
a full compensation of the resummation effects 
cannot be excluded. 

To draw firm conclusions on the small- X evo- 
lution of singlet structure functions also the next 
less singular terms have to be calculated. Since 
contributions even less singular than these ones 
may still cause relevant corrections, it appears 
to be indispensable to compare the correspond- 
ing results to those of future complete three-loop 
calculations. 
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